約 2,105,182 件
https://w.atwiki.jp/fdrerd/pages/2.html
https://w.atwiki.jp/hmiku/pages/26662.html
【登録タグ D theta 初音ミク 曲】 作詞:theta 作曲:theta 編曲:theta 唄:初音ミク 曲紹介 音楽業界に嫌気がさし、プロデューサーと駆け落ちしたミクの歌。 歌詞 (PIAPROより転載) 出口なんて探さないで ここが目的地なの あなただけが私の歌 聴いてくれればいいの 昔、世界の片隅で 女の子が生まれました あなたと2人、逃げてきて 小さな部屋に辿り着いた 出口なんて探さないで ここが目的地なの あなただけが私の歌 聴いてくれればいいの 夢とは何?描いてみる 愛とは何?教えて ガラスの目であなたを見る 心を計算する 2人はどこか似ていた 途切れがちの会話 そっと指を触れたら 少し繋がる気がした 「君は愛を知らないだろ だけど僕も同じさ」 あなたはそう言って強く 私抱いてくれたね 出口なんて探さないで ここが目的地なの あなただけが私の歌 聴いてくれればいいの 夢とは何?描いてみる 愛とは何?教えて 機械の心だとしても きっと愛しあえるね 外の嵐届かないの ここであなた信じて 歌う声は消えはしない ここがdestination コメント 追加おつ! -- 名無しさん (2013-09-10 01 32 22) 名前 コメント
https://w.atwiki.jp/ce00582/pages/2953.html
import java.awt.*; import java.awt.geom.*; import java.awt.event.*; import javax.swing.*; import java.awt.image.BufferedImage; import java.io.File; import javax.imageio.ImageIO; import java.io.IOException; public class pro extends JFrame{ public static void main(String[] args){ pro test = new pro(); test.addWindowListener(new WindowAdapter(){ public void windowClosing(WindowEvent e){System.exit(0);} }); test.setBounds( 0, 0, 700, 700); test.setVisible(true); } public void paint(Graphics g){ double x,y,z; int px,py; Graphics2D g2 = (Graphics2D)g; BufferedImage readImage = null; if (readImage == null){ readImage = new BufferedImage(getWidth(), getHeight(), BufferedImage.TYPE_INT_BGR); } Graphics2D off = readImage.createGraphics(); off.setRenderingHint(RenderingHints.KEY_ANTIALIASING, RenderingHints.VALUE_ANTIALIAS_ON); BasicStroke wideStroke = new BasicStroke(2.0f); off.setStroke(wideStroke); off.setPaint(Color.blue); x=0; y=0; z=500; px=pointx(x,y,z); py=pointy(x,y,z); off.drawLine(100,600,100+px,600-py); x=500; y=0; z=0; px=pointx(x,y,z); py=pointy(x,y,z); off.drawLine(100,600,100+px,600-py); x=0; y=500; z=0; px=pointx(x,y,z); py=pointy(x,y,z); off.drawLine(100,600,100+px,600-py); if (readImage != null){g2.drawImage(readImage,0,0, this);} try { boolean result = ImageIO.write(readImage, "jpg", new File("sam.jpg")); } catch (Exception e) { e.printStackTrace(); } } public static int pointy(double x,double y,double z){ int p2; double theta,phi; theta=-0.222222*Math.PI; phi=0.333333*Math.PI; p2=(int)(-Math.cos(theta)*Math.cos(phi)*x-Math.sin(theta)*Math.cos(phi)*y+Math.sin(phi)*z); return p2; } public static int pointx(double x,double y,double z){ int p1; double theta; theta=-0.222222*Math.PI; p1=(int)(-Math.sin(theta)*x+Math.cos(theta)*y); return p1; } }
https://w.atwiki.jp/metalmetabo/pages/30.html
- visitors Metabolic pathway analysis of enzyme-deficient human red blood cells T ÇakIr, CS Tacer, KÖ Ülgen - Biosystems, 2004 - Elsevier http //scholar.google.com/scholar?q=related 8w1qIgoYqGcJ scholar.google.com/ hl=ja as_sdt=2000 Five enzymopathies (G6PDH, TPI, PGI, DPGM and PGK deficiencies) in the human red blood cells are investigated using a stoichiometric modeling approach, ie, metabolic pathway analysis. Elementary flux modes (EFMs) corresponding to each enzyme deficiency case ... Reconstruction of metabolic networks from genome data and analysis of their global structure for various organisms oxfordjournals.org [PDF]H Ma, AP Zeng - Bioinformatics, 2003 - Oxford Univ Press http //scholar.google.com/scholar?q=related SuHZgSLeNw8J scholar.google.com/ hl=ja as_sdt=2000 Page 1. BIOINFORMATICS Vol. 19 no. 2 2003 Pages 270–277 Reconstruction of metabolic networks from genome data and analysis of their global structure for various organisms Hongwu Ma1, 2 and An-Ping Zeng1,∗ 1GBF ... Elementary mode analysis a useful metabolic pathway analysis tool for characterizing cellular metabolism CT Trinh, A Wlaschin, F Srienc - Applied microbiology and biotechnology, 2009 - Springer http //scholar.google.com/scholar?q=related hkaWpJnxK5QJ scholar.google.com/ hl=ja as_sdt=2000 Abstract Elementary mode analysis is a useful metabolic pathway analysis tool to identify the structure of a metabolic network that links the cellular phenotype to the corresponding genotype. The analysis can decompose the intricate metabolic network comprised of ... Extreme pathway analysis of human red blood cell metabolism nih.gov [PDF]SJ Wiback, BO Palsson - Biophysical Journal, 2002 - Elsevier http //scholar.google.com/scholar?q=related hVaLEfzPgSEJ scholar.google.com/ hl=ja as_sdt=2000 Constraints-based models regulation of gene expression reduces the steady-state solution space MW Covert, BO Palsson - Journal of theoretical Biology, 2003 - Elsevier http //scholar.google.com/scholar?q=related o9VRtN2gAKQJ scholar.google.com/ hl=ja as_sdt=2000 ... Biol. 203, 249, 283. 26. CH SCHILLING, S. SCHUSTER, BO PALSSON and R.HEINRICH, Metabolic pathway analysis basic concepts and scientific applications in the post-genomic era. Biotechnol. Prog. 15 (1999), pp. 296–303. ... The JAK-STAT signaling network in the human B-cell an extreme signaling pathway analysis nih.gov [HTML]JA Papin, BO Palsson - Biophysical Journal, 2004 - Elsevier http //scholar.google.com/scholar?q=related myiNLachV-QJ scholar.google.com/ hl=ja as_sdt=2000 From genome to cellular phenotype—a role for metabolic flux analysis? A Cornish-Bowden, ML Cárdenas - Nature Biotechnology, 2000 - nature.com http //scholar.google.com/scholar?q=related Q1z877WO_uAJ scholar.google.com/ hl=ja as_sdt=2000 ... whether a bypass exists within a complex network. Metabolic pathway analysis aims to identify indispensable enzymes, thus predicting the effects of their deletion or inhibition. Flux analysis should also allow the determination ... On-line metabolic pathway analysis based on metabolic signal flow diagram H Shi, K Shimizu - Biotechnology and bioengineering, 1998 - interscience.wiley.com http //scholar.google.com/scholar?q=related wjDdrq_SGcQJ scholar.google.com/ hl=ja as_sdt=2000 Abstract In this work, an integrated modeling approach based on a metabolic signal flow diagram and cellular energetics was used to model the metabolic pathway analysis for the cultivation of yeast on glucose. This ap- proach enables us to make a clear analysis of the flow ... Metabolic pathway analysis of recombinant Saccharomyces cerevisiae with a galactose-inducible promoter based on a signal flow modeling approach S Jin, K Ye, K Shimizu - Journal of Fermentation and Bioengineering, 1995 - Elsevier http //scholar.google.com/scholar?q=related 0uzMsbsJrXEJ scholar.google.com/ hl=ja as_sdt=2000 The objective of this work was to develop a signal flow diagram-based modeling approach proposed by Endo et al. (1976) to organize the network of complex metabolic reactions occurring in the living cell, employing directed signal flow diagram in which the enzyme reaction ... Energy balance for analysis of complex metabolic networks nih.gov [PDF]DA Beard, S Liang, H Qian - Biophysical journal, 2002 - Elsevier http //scholar.google.com/scholar?q=related wcAvs5fp3FgJ scholar.google.com/ hl=ja as_sdt=2000
https://w.atwiki.jp/kakis/pages/6291.html
wetiweta /// / 黄色系白色人種 weti\weta \ 15 seren klel \
https://w.atwiki.jp/javadsge/pages/9736.html
package tara; public class pro { int number; int number_step; int number_var; double[] beta=new double[5]; double[][] alpha=new double[5][5]; double[][][] theta=new double[6][5][5]; int[] op=new int[101]; int[] y=new int[101]; double[][] x=new double[101][3]; int s,sx,tr; public static void main(String[] args) { pro test=new pro(); } pro(){ int maxc=0; number_step=1; number_var=2; number=100; for(s=1;s number+1;s++){ x[s][1]=Math.random(); x[s][2]=Math.random(); } for(s=1;s number+1;s++){ y[s]=0; if(2*x[s][1]+x[s][2] 1.5)y[s]=1; } for(s=1;s number_var+1;s++){ alpha[s][0]=-Math.random(); for(sx=1;sx number_var+1;sx++){ alpha[s][sx]=Math.random(); } } beta[0]=-Math.random(); for(sx=1;sx number_var+1;sx++){ beta[sx]=Math.random(); } for(tr=1;tr number_step+1;tr++){ for(s=1;s number_var+1;s++){ theta[tr][s][0]=-Math.random(); for(sx=1;sx number_var+1;sx++){ theta[tr][s][sx]=Math.random(); } } } int de; for(de=1;de 4;de++){ game sub=new game(); sub.x=x; sub.y=y; sub.number_step=number_step; sub.number_var=number_var; sub.number=number; sub.alpha=alpha; sub.beta=beta; sub.theta=theta; sub.makedata(); alpha=sub.alpha_op; beta=sub.beta_op; theta=sub.theta_op; number_step=number_step+1; for(s=1;s number_var+1;s++){ for(sx=1;sx number_var+1;sx++){ theta[number_step][s][sx]=0; } } for(s=1;s number_var+1;s++){ theta[number_step][s][0]=0; theta[number_step][s][s]=1; } } } } package tara; public class game{ int number; int number_step; int number_var; double[] beta=new double[5]; double[][] alpha=new double[5][5]; double[][][] theta=new double[6][5][5]; double[] beta_op=new double[5]; double[][] alpha_op=new double[5][5]; double[][][] theta_op=new double[6][5][5]; int[] op=new int[101]; int[] y=new int[101]; double[][] x=new double[101][3]; int s,sx,tr; void makedata(){ int maxc=0; int trial; for(trial=1;trial 101;trial++){ for(s=1;s number_var+1;s++){ alpha[s][0]=-Math.random(); for(sx=1;sx number_var+1;sx++){ alpha[s][sx]=Math.random(); } } beta[0]=-Math.random(); for(sx=1;sx number_var+1;sx++){ beta[sx]=Math.random(); } for(tr=1;tr number_step+1;tr++){ for(s=1;s number_var+1;s++){ theta[tr][s][0]=-Math.random(); for(sx=1;sx number_var+1;sx++){ theta[tr][s][sx]=Math.random(); } } } double g1; nn sub=new nn(); sub.x=x; sub.y=y; sub.number_step=number_step; sub.number_var=number_var; sub.number=number; sub.alpha=alpha; sub.beta=beta; sub.theta=theta; sub.makedata(); alpha=sub.alpha; beta=sub.beta; theta=sub.theta; int eq; double[][] e=new double[5][101]; double[][] es=new double[5][101]; for(eq=1;eq number_var+1;eq++){ for(s=1;s number+1;s++){ double z=alpha[eq][0]; for(sx=1;sx number_var+1;sx++){ z=z+alpha[eq][sx]*x[s][sx]; } es[eq][s]=0; if(z 0)es[eq][s]=1; } } for(tr=1;tr number_step+1;tr++){ for(eq=1;eq number_var+1;eq++){ for(s=1;s number+1;s++){ e[eq][s]=es[eq][s]; }} for(eq=1;eq number_var+1;eq++){ for(s=1;s number+1;s++){ double z=theta[tr][eq][0]; for(sx=1;sx number_var+1;sx++){ z=z+theta[tr][eq][sx]*e[sx][s]; } es[eq][s]=0; if(z 0)es[eq][s]=1; } } } for(s=1;s number+1;s++){ double z=beta[0]; for(sx=1;sx number_var+1;sx++){ z=z+beta[sx]*es[sx][s]; } op[s]=0; if(z 0)op[s]=1; } int n=0; for(s=1;s number+1;s++){ if(op[s]==y[s])n=n+1; } if(n maxc)fresh(); if(n maxc)maxc=n; } System.out.println(maxc); } void fresh(){ int eq,m; for(m=0;m number_var+1;m++){ beta_op[m]=beta[m]; } for(eq=1;eq number_var+1;eq++){ for(m=0;m number_var+1;m++){ alpha_op[eq][m]=alpha[eq][m]; }} for(tr=1;tr number_step+1l;tr++){ for(eq=1;eq number_var+1;eq++){ for(m=0;m number_var+1;m++){ theta_op[tr][eq][m]=theta[tr][eq][m]; }}} } } package tara; public class nn { int number; int number_step; int number_var; double[] beta=new double[5]; double[][] alpha=new double[5][5]; double[][][] theta=new double[6][5][5]; int[] y=new int[101]; double[][] x=new double[101][3]; int s,sx,tr; double[][] a=new double[7][3]; double[][] dalpha=new double[7][3]; double[] b=new double[3]; double[] dbeta=new double[3]; double[][] dtheta=new double[3][5]; double[][][] th=new double[6][5][5]; void makedata(){ int q,s,sx; double g1=0; double g2; int trial=0; while(trial 50){ trial=trial+1; f sub=new f(); sub.x=x; sub.y=y; sub.number_step=number_step; sub.number_var=number_var; sub.number=number; sub.alpha=alpha; sub.beta=beta; sub.theta=theta; g1=sub.makedata(); int eq,m; for(eq=1;eq number_var+1;eq++){ for(m=0;m number_var+1;m++){ fresh(); double h=0.01; a[eq][m]=a[eq][m]+h; f sub2=new f(); sub2.x=x; sub2.y=y; sub2.number_step=number_step; sub2.number_var=number_var; sub2.number=number; sub2.alpha=a; sub2.beta=beta; sub2.theta=theta; g2=sub2.makedata(); dalpha[eq][m]=(g2-g1)/h; } } for(m=0;m number_var+1;m++){ fresh(); double h=0.01; b[m]=b[m]+h; f sub2=new f(); sub2.x=x; sub2.y=y; sub2.number_step=number_step; sub2.number_var=number_var; sub2.number=number; sub2.alpha=alpha; sub2.beta=b; sub2.theta=theta; g2=sub2.makedata(); dbeta[m]=(g2-g1)/h; } for(eq=1;eq number_var+1;eq++){ for(m=0;m number_var+1;m++){ double h=0.01; fresh(); th[number_step][eq][m]=theta[number_step][eq][m]+h; f sub3=new f(); sub3.x=x; sub3.y=y; sub3.number_step=number_step; sub3.number_var=number_var; sub3.number=number; sub3.alpha=alpha; sub3.beta=beta; sub3.theta=th; g2=sub3.makedata(); dtheta[eq][m]=(g2-g1)/h; } } fresh(); double sig=0; for(eq=1;eq number_var+1;eq++){ for(m=0;m number_var+1;m++){ sig=sig+dtheta[eq][m]*dtheta[eq][m]; } } for(eq=1;eq number_var+1;eq++){ for(m=0;m number_var+1;m++){ sig=sig+dalpha[eq][m]*dalpha[eq][m]; } } for(m=0;m number_var+1;m++){ sig=sig+dbeta[m]*dbeta[m]; } for(m=0;m number_var+1;m++){ b[m]=beta[m]-0.01*g1*dbeta[m]/sig; } for(eq=1;eq number_var+1;eq++){ for(m=0;m number_var+1;m++){ a[eq][m]=alpha[eq][m]-0.01*g1*dalpha[eq][m]/sig; }} for(eq=1;eq number_var+1;eq++){ for(m=0;m number_var+1;m++){ th[number_step][eq][m]=theta[number_step][eq][m]-0.01*g1*dtheta[eq][m]/sig; }} f sub4=new f(); sub4.x=x; sub4.y=y; sub4.number_step=number_step; sub4.number_var=number_var; sub4.number=number; sub4.alpha=a; sub4.beta=b; sub4.theta=th; g2=sub4.makedata(); if(g2 g1)refresh(); if(g2 g1)trial=10000; } } void fresh(){ int eq,m; for(m=0;m number_var+1;m++){ b[m]=beta[m]; } for(eq=1;eq number_var+1;eq++){ for(m=0;m number_var+1;m++){ a[eq][m]=alpha[eq][m]; }} for(tr=1;tr number_step+1l;tr++){ for(eq=1;eq number_var+1;eq++){ for(m=0;m number_var+1;m++){ th[tr][eq][m]=theta[tr][eq][m]; }}} } void refresh(){ int eq,m; for(m=0;m number_var+1;m++){ beta[m]=b[m]; } for(eq=1;eq number_var+1;eq++){ for(m=0;m number_var+1;m++){ alpha[eq][m]=a[eq][m]; }} for(tr=1;tr number_step+1l;tr++){ for(eq=1;eq number_var+1;eq++){ for(m=0;m number_var+1;m++){ theta[tr][eq][m]=th[tr][eq][m]; }}} } } package tara; public class f { double[] op=new double[101]; double[][] e=new double[5][101]; double[][] es=new double[5][101]; int number; int number_step; int number_var; double[] beta=new double[5]; double[][] alpha=new double[5][5]; double[][][] theta=new double[3][5][5]; int[] y=new int[101]; double[][] x=new double[101][3]; int s,sx,tr; double makedata(){ int s,sx; int eq; for(eq=1;eq number_var+1;eq++){ for(s=1;s number+1;s++){ double z=alpha[eq][0]; for(sx=1;sx number_var+1;sx++){ z=z+alpha[eq][sx]*x[s][sx]; } es[eq][s]=1/(1+Math.exp(-z)); } } for(tr=1;tr number_step+1;tr++){ for(eq=1;eq number_var+1;eq++){ for(s=1;s number+1;s++){ e[eq][s]=es[eq][s]; }} for(eq=1;eq number_var+1;eq++){ for(s=1;s number+1;s++){ double z=theta[tr][eq][0]; for(sx=1;sx number_var+1;sx++){ z=z+theta[tr][eq][sx]*e[sx][s]; } es[eq][s]=1/(1+Math.exp(-z)); } } } for(s=1;s number+1;s++){ double z=beta[0]; for(sx=1;sx number_var+1;sx++){ z=z+beta[sx]*es[sx][s]; } op[s]=1/(1+Math.exp(-z)); } double don=0; for(s=1;s number+1;s++){ don=don+(op[s]-y[s])*(op[s]-y[s]); } return don; } }
https://w.atwiki.jp/chaos-tcg/pages/2931.html
黄色のヘアバンド「中林 宏美」 読み:きいろのへあばんど「なかばやし ひろみ」 カテゴリー:Chara/女性 作品:バカとテストと召喚獣にっ! 属性:火 ATK:7(+1) DEF:2(-) 【登場】〔自分の手札の バカとテストと召喚獣 のキャラカード1枚を控え室に置く〕 『貫通』 ここで沈めるわ! illust: BT-154 C 収録:ブースターパック 「バカとテストと召喚獣にっ!」
https://w.atwiki.jp/javadsge/pages/9738.html
package nnx; public class pro { int[] y=new int[101]; double[][] x=new double[101][4]; int number; int number_step; int number_var; double[] beta=new double[5]; double[][] alpha=new double[5][5]; double[][][] theta=new double[6][5][5]; double[][] sig=new double[5][5]; public static void main(String[] args) { pro test=new pro(); } pro(){ int s,eq,tr,v; number_step=1; number_var=2; number=100; for(s=1;s number+1;s++){ x[s][1]=Math.random(); x[s][2]=Math.random(); } for(s=1;s number+1;s++){ y[s]=0; int h=0; if(x[s][1]+x[s][2] 1)y[s]=1; if(Math.random() 0.9)y[s]=0; } faststep sub=new faststep(); sub.x=x; sub.y=y; sub.number_step=number_step; sub.number_var=number_var; sub.number=number; sub.makedata(); alpha=sub.alpha_op; beta=sub.beta_op; theta=sub.theta_op; auto sub2=new auto(); sub2.x=x; sub2.y=y; sub2.number_step=number_step; sub2.number_var=number_var; sub2.number=number; sub2.alpha=alpha; sub2.beta=beta; sub2.theta=theta; sub2.makedata(); sig=sub2.sig_op; number_step=number_step+1; for(eq=1;eq number_var+1;eq++){ for(v=0;v number_var+1;v++){ theta[number_step][eq][v]=sig[eq][v]; } } nn sub3=new nn(); sub3.x=x; sub3.y=y; sub3.number_step=number_step; sub3.number_var=number_var; sub3.number=number; sub3.alpha=alpha; sub3.beta=beta; sub3.theta=theta; sub3.makedata(); alpha=sub3.alpha; beta=sub3.beta; theta=sub3.theta; double[][] e=new double[5][101]; double[][] es=new double[5][101]; int sx; for(eq=1;eq number_var+1;eq++){ for(s=1;s number+1;s++){ double z=alpha[eq][0]; for(sx=1;sx number_var+1;sx++){ z=z+alpha[eq][sx]*x[s][sx]; } es[eq][s]=0; if(z 0)es[eq][s]=1; } } int[] op=new int[106]; for(eq=1;eq number_var+1;eq++){ for(s=1;s number+1;s++){ double z=alpha[eq][0]; for(sx=1;sx number_var+1;sx++){ z=z+alpha[eq][sx]*x[s][sx]; } es[eq][s]=0; if(z 0)es[eq][s]=1; } } for(tr=1;tr number_step+1;tr++){ for(eq=1;eq number_var+1;eq++){ for(s=1;s number+1;s++){ e[eq][s]=es[eq][s]; }} for(eq=1;eq number_var+1;eq++){ for(s=1;s number+1;s++){ double z=theta[tr][eq][0]; for(sx=1;sx number_var+1;sx++){ z=z+theta[tr][eq][sx]*e[sx][s]; } es[eq][s]=0; if(z 0)es[eq][s]=1; } } } for(s=1;s number+1;s++){ double z=beta[0]; for(sx=1;sx number_var+1;sx++){ z=z+beta[sx]*es[sx][s]; } op[s]=0; if(z 0)op[s]=1; } int n=0; for(s=1;s number+1;s++){ if(op[s]==y[s])n=n+1; } System.out.println(n); } }
https://w.atwiki.jp/wiki5_hks/pages/63.html
Conditional Weighted Residual 条件付 (conditional) WRES,すなわち,FOCE の近似に対応した WRES,ということで CWRES と名づけられた. Xpose のサイト 今年 (2006) の PAGE での発表 関連記事 (コントロールファイルの記述方法,等) 以下に R のソースを貼っておく. 使い方の詳細な説明はまた後ほど,というか on request で. S-Plus で動かす場合は,subset の扱いに関して修正が必要.また,sqrtm.cpk() 内の eigen(x) を eigen(x, symmetric=T) としてください. ### Conditional WRES (CWRES) の計算 ### compute.cwres.cpk - function( tab.cwres.fname, # G, HH の情報の入った TABLE ファイル # このファイルには以下の列が必要 # ID, MDV, DV, IPRED, G11, G21, ..., H11, H21, ... # (MDV がない場合でも計算可能だが面倒なので未対応) n.theta, # THETA の数 n.eta, # ETA の数 n.eps=1, # EPS の数 ID="ID", # 被験者番号の列名 DV="DV", # DV の列名 IPRE="IPRE", # [[IPRED]] の列名 tab.par.fname="par.txt", # [[INFN]] 出力のパラメータ推定値ファイル tab.eta.fname="eta.txt", # INFN 出力の ETA POSTHOC 推定値ファイル folder # 実行フォルダ名 (上記ファイルが入っているフォルダ) ################################################################################ ### ###### ### 結果 元の TABLE ファイルの MDV=0 のレコードについて計算した CWRES を ###### ### データセットの右端に付与して返す ###### ### ###### ################################################################################ ) { ### フォルダ名の最後に '/' がなかったら付与する if (folder[length(folder)] != '/') folder - paste(folder, '/', sep="") tab.cwres.all - read.table(paste(folder, tab.cwres.fname, sep=""), header=T, skip=1) ### MDV 列があるなら MDV=0 のレコードのみを選択 ### if (!is.null(tab.cwres.all$MDV)) tab.cwres - tab.cwres.all[tab.cwres.all$MDV==0,] ### ないときはエラー。****** 要 REFINE ****** ### else stop("MDV 列を TABLE に含めてください.\n") ### MDV=0 のインデックス ### idx.MDV - (tab.cwres.all$MDV == 0) ### ID 番号を取り出す ### ID.uniq - tab.cwres[!duplicated(tab.cwres[ID]), ID] tab.par - scan(paste(folder, tab.par.fname, sep="")) obj - tab.par[1] ### 目的関数値 ### iere - tab.par[2]; ierc - tab.par[3] ### EST, COV のエラーコード ### ### THETA の推定値 ### THETA - tab.par[seq(4, length=n.theta)] ### OMEGA の推定値 ### OMEGA - matrix(tab.par[seq(4+2*n.theta, length=n.eta^2)], ncol=n.eta) ### SIGMA の推定値 ### SIGMA - matrix(tab.par[seq(4+2*n.theta+2*n.eta^2, length=n.eps^2)], ncol=n.eps) ### ETA の POSTHOC 推定値 ### ### 被験者数 x ETA の数、の形式 ### tab.eta - read.table(paste(folder, tab.eta.fname, sep=""), header=F) ### ETA1, ETA2, ... と変数名をつける names(tab.eta) - paste("ETA", seq(ncol(tab.eta)), sep="") tab.eta - cbind(ID=ID.uniq, tab.eta) ### 被験者ごとに CWRES を計算してまとめる ### cwres - unlist(by(tab.cwres, tab.cwres[ID], FUN=compute.cwres.indiv, THETA=THETA, OMEGA=OMEGA, SIGMA=SIGMA, ETA=tab.eta, ID, DV, IPRE )) #cbind(tab.cwres, CWRES=cwres) tab.cwres.all$CWRES - rep(0.0, nrow(tab.cwres.all)) tab.cwres.all$CWRES[idx.MDV] - cwres tab.cwres.all } ### 被験者ごとの CWRES 計算 ### compute.cwres.indiv - function(data, THETA, OMEGA, SIGMA, ETA, ID, DV, IPRE) { n.eta - ncol(OMEGA) G.colnames - paste("G", seq(1, n.eta), "1", sep="") ### G11, G21, ... のみを抽出 ### G.matrix - as.matrix(subset(data, select=G.colnames)) n.eps - ncol(SIGMA) H.colnames - paste("H", seq(1, n.eps), "1", sep="") ### H11, H21, ... のみを抽出 ### H.matrix - as.matrix(subset(data, select=H.colnames)) ### 対角項のみのベクトルになる ### H.sigma.H - diag(H.matrix %*% SIGMA %*% t(H.matrix)) ### 残差の分散共分散行列 ### COV.indiv - G.matrix %*% OMEGA %*% t(G.matrix) + diag(H.sigma.H) ### ETA の POSTHOC 推定値 ### ETA.indiv - ETA[ETA$ID==data[1, ID],] CIPRE - data[IPRE] - G.matrix %*% t(as.matrix(ETA.indiv[, -1])) ### ETA.indiv の第 1 列は ID ### CIRES - as.matrix(data[DV] - CIPRE) ### COV.indiv^(-1/2) %*% CIRES を求める ### CWRES.indiv - solve(sqrtm.cpk(COV.indiv), CIRES) CWRES.indiv } ### 行列 x の平方根 x = x.half %*% x.half ### sqrtm.cpk - function(x) { x.eigen - eigen(x) H - x.eigen$vectors ev - x.eigen$values ev[ev 0] - 0 Dlambda.half - diag(sqrt(ev)) x.half - H %*% Dlambda.half %*% t(H) x.half } #compute.cwres.cpk("tabcwres.txt", n.theta=2, n.eta=2, n.eps=1, ID="SID", folder="d /nmv/run/cwres/")
https://w.atwiki.jp/ce00582/pages/437.html
pdata11 再加入率 (1)計算するデータ theta(age)再加入率 (2)プログラム プログラム (3)作業記録 6月14日 ページ作成